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Abstract 

We derive the TBA system of equations from the S-matrix describing integrable massive 
perturbation of the coset Gi x G m /Gi +m by the field (1, l,adj) for all the infinite series of 
. simple Lie algebras G = A n , B n ,C n , D n . In the cases A n ,C n , where the full S-matrices are 

' known, the derivation is exact, while the B n ,D n cases dictate some natural assumption about 

the form of the crossing-unitarizing prefactor for any two fundamental representations of the 
algebras. In all the cases the derived systems are transformed to the corresponding functional 
^? F-systems and shown to have the correct high temperature (UV) asymptotic in the ground 

f--*) ' state, reproducing the correct central charge of the coset. Some specific particular cases of 

the considered S-matrices are discussed. 

^ 1 Introduction 

Thermodynamic Bethe Ansatz is know to be one of the most impressive achievements of two 
dimensional physics relating (in the high temperature limit) the data of a perturbed conformal 
field theory (CFT) and an exact factorizable S-matrix for relativistic integrable model with its 
spectrum and analytical structure. Unfortunately not only S-matrix often has a conjectural 
status, but also the TBA system corresponding to the S-matrix (e.g. [1] [2] [3] [4] [5]). The main 
obstacle in exact derivation of TBA system from S-matrix is usually related to the problem of 
diagonalization of transfer matrix (TM), which is especially a non trivial issue, when S-matrix 
has internal degrees of freedom, i.e. is non diagonal. In this paper we show how one can solve this 
technical problem for relativistic integrable models using some facts, established in investigation 
of trigonometric TMs of lattice models with Lie algebraic symmetries. These results were obtained 
in the framework of algebraic Bethe ansats [6] [7] , TM functional relations and analytical Bethe 
Ansatz [8]. 

The method to incorporate lattice model as a "carrier" of magnonic degrees of freedom re- 
sponsible for entire symmetries of a relativistic integrable model, is an alternative to another 
method of explicit lattice light cone regularization of whole relativistic model (see, e.g. [9] [10]). 
As far as we know, the first time this program was successfully brought about by Hollowood 
in [11], where he derived TBA system in the A n case. Another successfull implementation of 
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this procedure was done in [12] and [13]. In this paper we use this method applying it to the 
integrable models with other Lie algebraic symmetries. 

It is known for a long time [14] that integrable quantum field theories arising as perturbations 
of the coset CFT G Q* Gm by the operator (1, l,adj) for different Lie algebras G is a wide class of 

models, with with U q {G^>) symmetric trigonometric S-matrices , q a root of unity. More precisely, 
it was conjectured in [14] that the massive S-matrix has the symmetry U qi (G^) <g> U q2 (G^), 
where q\ = —e~ ln ^ l+9 \qi = —e~ ln ^ m+g \ where g is the dual Coxeter number of G. In other 
words, the S-matrix is, up to a scalar factor, a tensor product of two trigonometric U q invariant 
integrable models at q a root of unity, i.e. RSOS models. This class of S-matrices contains some 
other integrable quantum field theories (IQFT) with Yangian symmetry (rational S-matrices) , as 
a subclasses with specific choice of parameters l,m: Principal Chiral models (PCM): l,m — > oo, 
Gross-Neveu models (GN) 1 and current-current perturbations of WZW models: m = 1,1 — > oo. 
There is a lot of literature devoted to investigation of these models and we won't cite it here. 
In the context of the general U qi (G^) (8> U q2 (G^) symmetry these limits for m,l mean some 
essential simplification in the S-matrix structure, related to peculiarity of RSOS models and some 
identities existing between them for low level of restriction, which we discuss in the last section. 

As we said, the main goal of this work is direct derivation of TBA equations from the full 
S-matrix and their high temperature analysis. In principle, such derivation is impossible, in B n 
and D n cases, since the spectral decomposition of the S-matrices is unknown for arbitrary two 
fundamental highest weights of the algebras because of multiplicities appearing in irreducible rep- 
resentations decomposition. But as we will see below, even in these cases some information about 
the full S-matrix can be extracted from requirement of consistency of derived TBA equations. 
More precisely, the requirement that the obtained TBA system will correspond to the proper 
Y-system, fixes the crossing-unitarising prefactor of corresponding S-matrices. In contrast, in 
the cases A n and C n the irrep. decomposition of tensor product of any two fundamental repre- 
sentations is known, and full S-matrix may be written explicitly. Derivation of TBA Y-system is 
exact in the sense the assumption about this prefactor can be checked exactly. In all the cases 
of the considered Lie algebras G we show that the derived TBA equations lead to the correct 
ground state free energy reproducing in the UV limit the central charge of the coset. 

In the second section we show how results about transfer matrix diagonalization may be used 
in the derivation of TBA equations in the framework of Bethe Ansatz (BA) string hypothesis. We 
show the role of magnonic degrees of freedom and their relation to the main massive particles. In 
the third section we show that in the thermodynamic limit one of the possible magnonic degrees 
of freedom is always "frozen" - it always has zero density of holes, and we perform the reduction 
of this degree of freedom in the equations. The main fourth part is devoted to the transformation 
of obtained system of equations to the form of Y-system. As we show, this transformation is 
possible and requires a natural assumption about the form of the full S-matrix. Using results 
of [15], we show that the obtained Ysystems has correct high temperature behavior, reproducing 
central charges of the coset for any G = A n , B n ,C n , D n . In the fifth section we show that the 
assumption made about the form of the crossing unitarising scalar prefactor of the S-matrix, is 
really correct in the A n C n cases. We also discuss particular cases of the general S-matrix for 
some specific models with Yangian symmetry. We conclude the paper by brief discussion. In 
Appendix we collect kernels and technical details of the TBA derivation. 

1 Note that not any Lie algebraic symmetry may be a symmetry of GN model: its impossible to construct a 
GN-like interaction of Majorana fermions with Sp group. 
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2 Transfer matrix diagonalization and TBA equations 



In contrast to the case of relativistic two dimensional integrable quantum field theories (IQFT) 
with elastic (diagonal) S-matrices, such as, e.g. S-matrices of affine Toda field theories (see, 
e.g. [16] [17] [18]), the transfer matrix diagonalization for the IQFT with internal Lie algebraic 
symmetry, is hard, and, in general, non solvable problem. More progress was achieved in solution 
of this problem in lattice models, such as spin chains or RSOS models invariant under some Lie 
algebraic symmetry, rather in relativistic IQFT. Here we recall one simple method to attach the 
lattice results, derived in RSOS like traditional Bethe ansatz methods, to explicit derivation of 
TBA equations for S-matrices describing integrable perturbations of G ' xGm coset CFTs (G = 
A n , B n ,C n , D n ) relevantly perturbed by operator (1,1, adj). The S-matrix for this IQFT was 
conjectured a long time ago [14] and it has the form 

s ab (e) = x ab (e)sS(e)®s^\e), (i) 

where S^(9) is unitary, crossing symmetric "minimal" (without poles in the physical strip < 
Im9 < 7r) RSOS-like S-matrix for scattering of two particles from two multiplets corresponding 
to fundamental weights a and b of algebra G. k is the restriction level of RSOS model. X a b is a 
CDD factor which generates poles for the S-matrix corresponding to each fundamental weight of 
G, and guarantee the bootstrap closure. Recall that from the point of view of quantum groups, 
Sq)(9) RSOS S-matrix has the U qi (G^) symmetry with q a root of unity q = —e~ l7T /( l+9 \ 

The procedure of TBA derivation is standard: we pull one particle j from the fundamental 
multiplet cij with rapidity 9j, through a gas of other particles living on a circle of length C. On 
the way it scatters on each other particle with rapidity 9i with the S-matrix S ajai (9j — 9i), 
giving rise to the transfer matrix 

M 

T*i{e j \0 il ,...,O iN )= II S ajai (9 j -9 i ). (2) 

The requirement of the wave function periodicity looks like 

^ !:,h "V" !^^., U X .U- 9 j ,; = 1. (3) 

Non diagonality of the scattering leads to the change of states of the particle inside its multiplet. 
In terms of Bethe ansatz, this change is taken into account by means of magnonic excitations 
described by Bethe ansatz equation (BAE). They are responsible for the non diagonal part of 
the S-matrix, defined by spectral decomposition, whereas the diagonal part is defined by the 
prefactors before this spectral decomposition - X a h, and crossing- unitarising prefactors cr^ m ^ of 
the minimal S-matrices S^ m \ Explicit and full form of the non diagonal part of the transfer 
matrix in terms of magnonic degrees of freedom is complicated. But it was proven by Kirillov 
and Reshetikhin for simply laced algebras [6], and conjectured, and partly proved, for non simply 
laced algebras [8], that in the thermodynamic limit, when the number of particles Af together with 
the length C are going to infinity, there is dominating "top" term for transfer matrix. Leaving 
only this top term, we have for the transfer matrix the following expression 

M 

T^(9 j \9 il ,...,6 iN )= [J X ajai (9 j -9 i )a^ ai (9 j -9 i )a^ i (9 j -9 i )x (4) 
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x ^ Slnh (TO - g + <)) ^ Sinh (gngg (gj - g + ''a/)) 

«=i sinh ( i " n ° " ^/)) ° =1 sinh (sp^y " v * ~ 

Here the last line describes the "top" term contribution, according to [8], g is a dual Coxeter 
number of algebra G, t a - integer number related to the node a of Dynkin diagram of G. For 
algebras G = A n , B n , C n , D n considered in this paper, it is equal to 1 for long root nodes, and 2 
- for short root nodes. Sets of numbers Ua , Va satisfy the BAE [19] 



sinh {wkj K - °j + iu}a > - a b)) _ nb » ^ sinh (wk) (j - g + • a 0) (5) 

i=i sinh (lii+gj ( u a ~ 6 i - ■ a b)) c=l (3=1 sinh {^{jTg) (< - Up - i(*b ■ Oi c JJ 

and the same for v^, with I replaced by m. Here uj a ,a a are fundamental weights and simple roots 
of G. is a constant which is not important for us here. 

It is worth to notice here that rapidities of the physical particles 6j appear as inhomogeneities 
in the l.h.s. of the BAE © for magnonic degrees of freedom. This procedure differes from light 
cone lattice regularization scheme for relativistic IQFT, when one considers BAE like (J5J) with 
light cone inhomogeneities G in the l.h.s., and mass scale is introduced in a special scaling limit: 
G — > oo, lattice step a — > 0. 

It is important that according to the general conjecture [6] [8] about the top term (J3J, this 
transfer matrix eigenvalue is associated to a representation of G n x G n with highest weight 

( M n M n \ 

= £^-E M iV>E^ -£MMa , (6) 

\i=l a=l i=l a=l / 

and the RSOS restrictions n®6 < 1,^6 < m impose important restrictions on the possible 
values of Ma , where 9 is the highest root of G. 

Procedure of taking the thermodynamic limit j\f — > oo, C — > oo is standard: rapidities 6j 
become dense, and solutions of BAE form strings. The string hypothesis for any algebra G 
was formulated in [7] and it looks as follows. In the thermodynamic limit macroscopically large 
amount of solutions of BAE have the form of color a, fc-strings 

< = u a , k + it-\k + 1 - 2j), j = 1, .., k, (7) 

where u a ^ is real and has some density (density of string), and possible values of k are k = 
1, t a l. 

Before we start the thermodynamic calculation let us fix some notations. We use the following 
Cartan matrices for G: C a b = " a ^ h ■ As we already defined, t a = - 2 a . The fundamental weights 
satisfy a a oJb = ^ab^a 1 - We also define t a b = max(t a ,4). We use the following Fourier transform 
convention 

-I /'OO /'OO 

/(«) = 7T / /He^dw, f(u) = / f{u)e~^du. 

^ J —oo J —oo 

Now the standard thermodynamic calculation goes as follows (see, for example, [6]). Consider 
the logarithm of © . The possible ambiguity 2niq can have a holes in its continuous occupation if 
integer numbers q. In the thermodynamic limit such holes form hole densities. We sum up these 
equations over belonging to a color a fc-string, introducing densities pi for real coordinates of 
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strings ii a ,fc, and hole densities pf, for holes in 1/ C normalized q distributions. The same procedure 
we perform with the BAE for v^, introducing magnonic densities rji for real coordinates of strings 
v a ,ki and densities of holes rfe. The same procedure can be applied to the In of the eq. (jHJ) with 
explicit form of the transfer matrix given by (jlj , with introduction of particle a densities a a , and 
hole densities for them a a . The resulting equations have the form [7], which is the simplest after 
the passing to the Fourier transform. Here and below we will work mostly in the u space, so we 
will omit the hat on all the variables depending on u>, and moreover, will omit their argument co, 
except for the cases when it will be different from u. 



n t a l t a m 

JSfl - E W - E 4a! Vj - E «,%*■$ «» 

6=1 j=l j=l 

n t^l 

4tf v = #+£E^£* w V* ( 9 ) 

6=1 k=l 

n t^l 

a$*»o* = ^ + £E M «»4r 9)i U (10) 

6=1 fc=l 

Here a = 1, ...,n, j = 1, ...,t a l in eq. @, j = 1, ...,t a m in eq. I)1U[). Masses m a of multiplets for 
different algebras G will be listed later, cosh# is Fourier transform of cosh (9 function, and the 
kernels in the u space look like (here and in what follows we use the short notations for sinh and 
cosh functions sinh(wx) = [x], cosh(u;a;) = (x)) 

~ (11) 



A {L)jk 



[L] ' 

[mm(jt-\kt^)][L - maxjjt- 1 ,^ 1 )] 

Kim 



M ab = C ab ^ + 25 ab ((t- l )-l), 

tab 

and Y ab is Fourier transform of Y ab (9) which comes from the S-matrix S ab (0) prefactors 

The equations (|5|)- (|10|) are the basic equations for the TBA derivation. Before we come to it, 
one important step is necessary. As it was firstly noted by Bazhanov and Reshetikhin, one of the 
degrees of freedom in these equations is frozen. As we will see in the next section, oppositely to 
spin models, the RSOS restriction dictates jo" ; = rf® m = in the thermodynamic limit, which 
can be used for reduction of the highest t a l and t a m strings. 



3 Maximal string reduction 

Consider the zero mode of the /-system Q for j = t a l. Using explicit form of the kernels pijl. 
one has 

r-° a W = KM + j*- £ c ab J2 kp b k (oy 

y y 6=1 k=l 
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In the thermodynamic limit the highest weight (JBJ), which dominates in the transfer matrix 
eigenvalue, becomes 

n / t a i \ 

a=l \ k=l / 

and due to the previous equation, using Ylb=l Cab^b = ct a , it gives 

; _i_ n 

^ = ^E«°k- 

q 

y a=l 

RSOS restriction pfl)0 < I now gives ^yELiPyCM ^ L The fact that ma y be 

only non negative, and that uj a 9 is a positive number for any a and any G, necessarily requires 
that p1 a i(0) = for each a in C — > oo limit. From this immediately follows that p^j{uj) = for 
any cj. The same is valid for rj+ AoS) = 0. 

Using this fact, we express p^j through other variables. Eq. Q gives for j = t a l 

„(*+$) _a _ V* A/f 4 (I+p)t»I,fc J. , M \9\ M a6 6 

6=1 fc=l [ + 9i 6=1 raft J 

The inverse of the matrix , M jf , we will denote A ab n . We will use also another matrix A ab " = 
2W^G„_ Tlie inverse Q f foe matrices A G b n will be called K G b n = (A Gn )~* . A Gn can be calculated 

case by case for each of the four algebras we consider here. The list of matrices A G b for each G 
one can find in the Appendix. So we have 

1 n t a l—l 

pw=mE^ 6 -E41^ (13) 

L J 1 — 1 1 — 1 



In the same way one gets 



6=1 k=l 



ta.m—1 



c = R E^- E 41^- (14) 

1 ' 6=1 fc=l 

Substitution of (|13|) and ((Til) into © and ifTUj) gives, after some simple algebra, reduced magnonic 
BAE: 

a %y = ^+EE M ^if>£ (15) 

6=1 k=l 



n %l— 1 
6=1 k=l 



Before substitution of (|13|) and (|14|) into the massive equation (JSJ) we will make an important 
assumption: for all the algebras G the kernel Y ab has the form 

Y ab = i P A^, (17) 
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with some scalar function (p(u>). As we will see later, this assumption is correct in both G = A n 
and G = C n cases, when the full S-matrix, including its spectral decomposition, is known. With 
this assumption substitution gives 



t a l—l t a m—l 

,(0 „a J m ) 



cosm- Agy - E <U - E <W 



where 



^"^ + a ' [/]2(l) +flm [m]2(l)' J 



So, effective reduction of the thermodynamically "frozen" degrees of freedom leads to the system 
of equations (|18j l .(|15 |) and (|16j). very similar to the original ones. The effect of reduction is the 
change of parameters I + g —* I, (m + g — ► m) and possible string length is now not grater then 
t a l - 1, (t a m - 1). 



4 Thermodynamic and transformation to Y~-system 

In the massive BA equation (|18|) we have magnonic degrees of freedom - magnonic densities, 
which represent not excitations but rather Dirac vacuum for magnonic degrees of freedom. If 
we want the massive equation to be explicitly dependent on magnonic excitations, represented 
by hole densities Pj,rjj> we should solve eqs. (|T5|l and ifTfij) with respect to Pj,r]j as functions of 
Pj,Vj an d cr a , and substitute these solutions into the ea.qi8j). As we will see on the stage of doing 
thermodynamics, this will give us a possibility to convert the system ()18)) . (jl5) ). (jlf))l more easily 
to the form close to a so called Y-system. We have it as a goal, since such kind of systems were 
classified [15] according to affine symmetry standing behind the system, and this classification 
adjust to each such system its UV limit (in thermodynamic terms, the limit T — > oo ). This 

G x G 

limit contains information about the central charge of the unperturbed CFT - the coset —k — — . 
Actually, we don't need the expressions for Pj,r]j themselves, but only the specific combinations 

entering the massive equation, like Y^=i X a j/t Pv 



4.1 Transformation of magnonic degrees of freedom and S-matrix prefactor 

General description of the procedure is possible, but we will do calculation separately for both 
of non simply laced algebras, since this will make the formulas more transparent. We multiply 
eq. 1)15)1 by the matrix 

K k J = S kj + ^ry (C^ 1 - 1 - 25 kj ) = K^{t~% (20) 
and sum over j from 1 to t a l — 1. Using the identity 

tal-l 



one gets the equation 



y Kg? aft). = -Mr-, (21) 

^ a 3^ 2{ta 1 ) ^ ' 



2( - ta ) y-l h ! J 
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After some algebra the kernel J may be written in the universal form: 

ab 2(ta 1 )[t- b 1 ] Ua J tbKUj y ' 

t b /t a -i 

+ E ^6 ^ (^i(*+l)-*o9,*oJ + K(k-l)+t a q,t a j) 
q=l 

Simply laced cases (A n ,D n ). 
In these cases all t a = 1, and J has a simple form 



T kj _ M ab x 
J ab ~ 2(1)°^' 

and eq. (122H has the form 



c i— 1 n 

K ' i=i 6=1 

Multiplying it by aj? with summation X^i=i> an d using the identity (|21jl. one gets the equation 
for variable :c a = £]fc=i Pk' a = h--i n 



6=1 

which one can easily solve using the inverse kernels A Gn (|77j) or (|80j): 

1 n 

^ = 2^1^(4^-^). (25) 

^ ' 6=1 

Non simply laced cases. 

Unfortunately, the calculations in non simply laced cases are technically more involved, al- 
though straigtforward. We decided to present them in details. When some of t a are equal to 2, 
one has from (J23|) the following form of J: 



when t a = tb = 1, 

T kj _ M ab 
2(1 



J ab ~ n/n^J' 



when t a = 1, t& = 2, 



J. 



Ay 



(i) — 2 (i) 

when i a = 2, 4 = 1, 



^"2(1)^' 
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when t a = 2, % = 2, 



Case C n . 

Equations (|22j) can be written as 



r 21-1 „ n-l 

On 



2(1 V = E ^ + 2(172) g M -"? + 2072) M «"^ 

Here in the first equation a = 1, ...,n — 1, j = 1, 2/ — 1, and in the second - j = 1, Z — 1, and 
the kernel 

1 



G ' = ^ + 2(1/2) < 26 ) 

In these equations and below a fractional index at any variable p, p, rj, rj means that it is equal to 
zero. We would like to rewrite these equations using K kernels instead of M (see Appendix): 

= E^ fc (I)^+7E^>S + 7^ 2 (27) 

r. Z-l n-12l-l 

^-a* = Y,K ik (")r k +J2T, GjkK %p b k+p]- (28) 

^ ' k=l 6=1 fe=l 

Multiplying the first equation by ajj 2 and taking the sum Y^j=\ ■> an d the second - by 
and taking the sum we use the identity 

fri 3 k/2 2(1/2) 

and explicit form of M. We obtain the following system of equations 

n-l 



2(l/2)^/ 2 " 2(1/2) Pl + ^ af> +A « X 

2(1) Gl " 2(1) Pl (1) X +X + 2(l) Pl m 

for the variables x a = Y^j=i a f/2Pp a = l'-- n ~~ 1> an d ^ = X^=i a jVj- Their solution (for 
details see Appendix) has the form 

n-l 

c„ /JO 6 ~b 



(29) 



^(afV-^)+^(2(l/2)^( 



2(1) 

-i4& [a b -2(1/2) ^i^/2)^ 



a 1/2 a 



6=1 

21-1 



k=l 



(31) 
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r a _ (1/2) ( f {l) \ ^ / (I) c b Cn ^ 

X ~ (1) I 2(1/2) V 1 Pl J + ^ l a i/2^ a A* ft 

(^) _1 iS {^f/ - ^E^l fc (a;/2)p|j . (32) 

Case B n . 

In a similar way one can deal with the B n case. Eqs. (1221) in this case take the form 
-J^a a = Y, Kik WR + ^E M ^' a=l,...n-2 



2(1) ^ w ™ 2(1)^ 

w fe=l w 6=1 

^y^' 1 = E ^ H ^ + 271) K-l,n-2Pr 2 + Mn-Ln-lPr 1 ) + 



fc=l 

2Z-1 



f 2i-l 1 1 

2<tj2f n = £**<»/*>« + _M„,_ 1P ; 72 ' + — AW», 

where j = 1, ...Z — 1 in the first two equations, and j = 1, ...2/ — 1 - in the last one. G is the same 
as in eq. (|26jh In terms of kernels (see Appendix), these equations look like 



2(1) 



Z-l n 

E^'M^ + E^L a = l,...n-2 (33) 
j=i 6=1 

r 2-1 n-1 21— 1 

'»,»-*■ = E^'H^+E^-^+^E^i ( 34 ) 



fc=l 6=1 fe=l 

22-1 



2(1) 

Jhl_a n = E^'(^/2)^ + y<"n-iP^ 2 1 +^- (35) 

Multiplying (j"*"*"|) . ((3*1)1 by a^} with summation X*fc=i> an< ^ _ by a^j 2 with summation Y^k=\i 
as in the C n case, we get the following system of linear equations 



2(1) 2(1) 



1 n 

^TPi + E^ 6 ' o = l,...n-2 (36) 



6=1 

,(0 1 n ~ l 



'1 



2(1) - 2(1) ^ + g + A " --'-" " 2m K ^ (37) 

(0 

a 1 1 

1/2 n _ 1 -xn , 1 t^B n ,^n — 1 1 



2(1/2)'" = 2(172) « + +X " (38) 

with respect to the indeterminates x a = Y?j=i a j <P < ji a = I, —n — 1, and x n = Y^f=i a j/2Pj- 
Their solution (see Appendix) is 

1 1 A ab / (0 J ~b\ j_ (1/2) , Bn / (J) n ~„\ (1/2) 



£ ^ y^l " + (4/V n " P?) " ^fnVl, (39) 
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1 A nb (Al)h ~b\ , (V 2 ) A B n (SI) „n ~n\ A nn-1 n 



6=1 



E ^ (< v - ) + (4/V n - p?) - i^V • (40) 



The same solutions y a , y n one obtains for the other (rj) magnonic contributions to the massive 
equation (|18|) , They have the same form with change p — > rj, p — > rj, I — > m 

Now we substitute the obtained expressions for magnonic transfer matrix contributions (|2l))) . (f52*|) . l)31j) . p§|) . ()4Uj) 
into the massive equations (|TH|) . The non simply laced cases will again be considered separately, 
but we start from the simply laced cases. 

Simply laced cases (A n ,D n ). 

Substitution of ()25|) . and the similar expression for the rj dependent part, into ()18|) gives 

s- = ^55T* - + 2S^) | ^ + ^ | A- + a) . 

Multiplication of this equation by the matrix K Gn , the inverse one to , gives the equation 



6=1 



The main feature of the mass spectrums m a for all the algebras, listed in the Appendix, is that 
the vector ^-cosh# is an eigenvector of Kj^ with zero eigenvalue, which leads to disappearing 
of this term in the last equation. It makes possible to transform the system to so called universal 
form, when it does not contain any other external functions or parameters, but only the variables 
themselves (see below). 

Together with (|19j) the coefficient before a a in the last equation may be written as 

^ + 2(1) {[l) ai +Ql + [m] am +Ql 

which after some algebra can be written as 

Jl+9) j_ n (rn+g) 
I "1 ' 1 

— • 

The main conjecture is that 

tn — ( A A 2B+l+m-l\~ _ ( ..^g+i+m-A ~ f 

Y — y A g+ l t g + l J — \^g+ m ,g+m J • l 4Z J 

As we will see in the next section, this is exactly what one has from the S-matrix in the A n case, 
and we suppose the same expression is correct for D n case too. This assumption will be shown 
correct for S^ 1 . With this expression for ip, after some algebra, one gets the coefficient before 
a a equal to 1. 

Non simply laced cases. 

As we will see now, this property of corresponding coefficient will remain valid in non simply 
laced cases too, and will define the form of S-matrix prefactor ip. 

CclSG. 
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Equations (|18|) written separately for a < n — 1 and a = n, look like 

n 

a a = ^^Td-5^A c Ja b -x a -y a , a<n-l, 

6=1 

71 6=1 

Substitution of (|3*T|) and (|3*2*)l into them, after some algebra, using Y^l=i ^ab^-bc = ^ ac ' gi ves 
the following equations 



n 
6=1 



6=1 



(1 /|l(aW +a^)a»- (J/2) '^ 
2(1) + 2(1) 



VxJ-S* = _~ a n_ «i +ai ff . + ft+i + 
Z^i nb r 9 1 ^ 1 



,, , As V 2 d) (!) ,, 

+{p —> Tf,l—>m}. 

We recall that massive terms disappear since they are eigenvectors of K Cn with zero eigen- 
value. Using the ea.(|27|) at j = 1 and explicit form of the kernels, one can see that the last sum 
in the last equation is nothing but ^hjPi ~ > which gives the system 

ElfCn~b ( ~ , (V 2 ) / (0 , (rn)\\ a , (V 2 ) /-a , ~aN /, a \ 

ab = ~[<P + -^y- ^i/2 + a i/ 2 J J a + -^y- (Pi + »h ) (43) 

Vi^"? 6 - - 15+ 1 + 1 U n + Pl +Vl + ^ +7?1 (44) 

2_^ K nb° " 1^+ 2(1) I + 2(1) + 2(1) 



6=1 



We make here the same conjecture as in the simply laced cases, which will be checked by explicit 
calculation from the S-matrix in the next section 

V ~ { A 9+l,9+l 

(Recall that for C n g = n + 1). As we saw above, this form of ip leads to the fact that the 
coefficient before a n in eq. Q44|) is 1. Using this fact, one can easily see that the coefficient before 

2{ , l/2) 2 

a a in eq. (|43|) is equal to / = ^ fa ' . So the final form of equations in the C n case is 

-V ^K Cn a b - -a a + ^ + ^ (45) 
fh ab + 2(1/2) (45j 



6=1 



^ nb 2(1) 2(1) 



B n case. 
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Doing the same calculation as in the C n case, using l)3"5)l (110)) and the identity 

1 ' /^\ n,n— 1 J n tmi 



one gets 



l^K ab a = -\<p+ 2(1) Icr +__ + *„,„_!-_ (47) 



6=1 



where in the first equation a = 1, ...,n — 1. 

The conjecture about the form of </? 1)42)1 remains unchanged, which gives, as we saw, coefficient 
1 before a a in (|47j). and coefficient / before <r n in ()48|). It gives the following final form of equations 
in the B n case 

~ f l^K nb a = -a + Wm - (50) 
J 6=1 w ' 

We see that the equations we got has a compact form, and are similar in simply laced and 
non simply laced cases. Their universality is in particular expressed by the fact that they don't 
contain mass terms. 

Before we will do thermodynamic of the system we prefer to rewrite magnonic and massive 
equations as one equation. It can be done if we introduce the following new notations 



b 3 


= 


b 3 


= P% j = 1, 


it a l 1, 


3 a . 
b J 


= % 


■ S J 


= j = 'I, 


...,-t a l + l, 


• s o 


= cr , 


s 


= a a . 





As one can see, in terms of variables s, s, the role of holes and pseudoparticles flipped for magnonic 
degrees of freedom and remained the same for massive ones. One can see that in terms of s,s, 
simply laced equations 1)41)) and (|24j) together with the copy of (|24|) for 77, rj, can be written as 
one equation 

n l-l 

E K a b ^ = - E **W> (51) 
6=1 j=-m+l 

where G is either A or D, and the kernel K k i has the same definition as before (|2U|). but its 
indices are running from —m + 1 to I — 1. The equation with j = is now the massive equation. 
In the C n case equations 1)27)1 . its analog for 77,77 and (|43j). can be written as 

n-1 21-1 
7E^ + 7^/2 = - E ^(«V 2 K> a < n-l (52) 

6=1 ^ fc=-2m+l 
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Equations (|28[). their double for 77,77, and l|44|) one can write as the following one equation 

^"Tlf E C^-^- £ K jk (u) si. (53) 

^ ' fc=-2m+l fc=— ra+1 

In the same way the B n case equations (|49j) and (|33|) . (|34|) with their 77,77 analogs give 

j;^-^ J2 G^ = -Y d K kj {u)s a j ,a<n-l (54) 

6=1 ^ ' j=-2m+l j=l 

and equations (|5Uj) . (|35|) with 77 partners 

21-1 

^-Jm)^ = - E K^( U /2)sl (55) 

^ ' ' fc=-2m+l 

Doing thermodynamic is standard procedure (see for example [20], [21], [22]): we should 
minimize the free energy F = E — TS, (E - is energy, T - temperature, S - entropy), using the 
derived equations as constraints. The energy is the energy of massive particles 



E = ^2 d9s%(9)m a cosh 0, 

a=l J-oo 



and the entropy can be calculated from combinatoric of states as particles and holes in the 
thermodynamic limit 

71 t a l 1 /»00 

0=1 j=-t a m+l J -°° 

We will not repeat the standard free energy minimization procedure. If the starting constraints 
equations were written in the form 

n t a l—l 

E^ S 5 = - E pjks ^ 

6=1 k=— t a m+l 

with some kernels Q, P, as it was in eqs. (|51 |1 -(|55 |) . the variation leads to the set of equations in 
the form so called y-system 

n t a l—l 

E^W+b = E p3kL U^ ( 56 ) 

6=1 k=— tam+l 

where we defined "dressed energies" e|, and 

4 = exp(4) = Y k a , L a (±)k = In (1 + exp(±4)) . 
s k 

The free energy itself can be expressed in terms of the stationary values of functions Y k a (±oo) in 
the rapidity space. The well known relation between the central charge of the relativistic theory 
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defined on a cylinder in the UV limit (radius of the cylinder goes to zero, when temperature is 
going to infinity) F = — ^T, gives a possibility to extract the central charge of the corresponding 
CFT in the UV limit ( G ' xGm in our case) 

7T 2 ^ 

a 

Here L is the dilogarithm function, Yj,Y™ are the stationary values of corresponding variables, 
introduced above, for the full system, and for the system with removed massive variables Yq. 
This result can be considered as the most serious check of the S-matrix, and the way from S- 
matrix to the central charge usually has status of conjecture in the literature, dealing with TBA 
analysis. 

Explicit use of the kernels form for the systems l|51j l — Q55 |) after their transformation into the 
thermodynamic ones according to (|56|). gives the following y-systems in the rapidity space. 





(57) 



D, 



r^o) = (1 + (1 + Yf-^ef) 



R%(6) = (l + >7 +1 (0)) (l + ^(A)) , a < n - 3, 

R]- 2 (6) = (l + Y™- 3 (#)) (l + F/- 1 ^)) (1 + Y?(9)) , 

R^-\9) = l + Yf- 2 (9), 

iqip) = \ + Y^- 2 {e). 



*m = (l+^L^iW)(l+ir X W),a<n-l, 

R]{9) = (l + Y 2 y\9 + i/2)) (l + Y 2 y\9 - i/2)) (l + Y 2 n r + \(9)) (l + Y%\{jf)) . 



r*(P) = (i + y/ +1 (0)) (l + Y*-\9)) * 

[(1 + Y%{9 + i/2)) (1 + Y%{9 - i/2)) (1 + Y 2 n j+1 (9)) (l + ^(fl))]'*' 
R«(9) = (l + Yf /2 \9)). 
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Figure 1: TBA diagrams for the cases a (left) G = A 5 ,l = 4, m = 3, b (right) G = B 5 ,l 
3, m = 2 
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Figure 2: TBA diagrams for the cases a (left) G = C§,1 = 2, m = 2, b (right) G = D-j,l 
3,m = 2 



Here 



1 + 



i + 



+ i/t a )K a (0-V*a). 



One can see that left hand side Rj of these equations contains different j indexes and the same 
index a, whilst their right hand side contains different a indexes. Traditionally these TBA 
equations one represents schematically as TBA diagram (see fig. 1,2). Their nodes correspond 
to each Y?. Those j's which appear in the Rj along with j itself, are depicted by vertical lines 
connecting (a, j) node with others. Horizontal and other lines represent other indices then a and 
j appearing in the r.h.s. of (a,j) equation. On these figures massive nodes are depicted as bold 
ones, and magnonic - as grey ones. The purely magnonic equations are represented by the same 
TBA diagrams with removed massive nodes. 

As one can check, the above Y-systems exactly reproduce those which were considered in [15], 



if one changes Y? 



3 



In [15] such y-systems were related to quantum affine Lie algebras 



U g (Gn ) and their representations. The mathematical meaning standing behind functional Y- 
systems and their relation to representation theory of quantum affine Lie algebras, their character 
relations and transfer matrix functional relations, is deep and beautiful issue, requiring further 
investigation. Here we will cite the main result of [15]: the dilogarithm sum rules like (|57|) 
corresponding to the ground state of the y-systems we obtained (we omit some technical details 
and present it in a simplified form relevant for our case): 




/ dim G T: 
1 + 9 



n, 



(58) 
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where X® satisfy the same functional relations as our [Yj'j ■ Using the relation L(t) = 1 — 

L(l — t) and substituting into 1)57(1 . one can see that the latter reproduce the correct central 
charge 

„ ( I m Z + m \ 

c = dim G n 1 

\l + g m + g l + m + gj 

for all the coset models considered here. 

This completes the derivation of TBA y-systems from the S-matrices of relativistic integrable 
models, and their check by comparison of their ground state thermodynamic with the central 
charge of their UV limit CFTs. This is one of the most convenient checks of the S-matrix 
correctness. 



5 Full S-matrix 



As we have seen in the previous section, the consistent TBA, which reproduces the correct central 
charge, dictates unique form of the S-matrix crossing-unitarizing prefactor together with CDD 
facor, for any pair of two fundamental representations for any algebra G. Here we will show that 
in the cases where the full S-matrix is known for any pair of fundamental representations, i.e. 
in multiplicity free cases A n ,C n , this prefactor can be calculated from the S-matrix. In the A n 
case this calculation was done in [11]. The most interesting information one can extract from the 
above TBA analysis, is the prefactor for B n ,D n cases. Although the S-matrix S ab is not know 
in this cases for any a, b, it is known, for example, for spinor-spinor S-matrices: in B n case for 
a = b = n, or in D n case - a = b = n and a = n, b = n — 1, since then spectral decomposition 
is multiplicity free. Corresponding rational S-matrices were constructed for PCM in the seminal 
work [23], but we, unfortunately, don't know any published trigonometric generalization of their 
result. It would be interesting to work out corresponding trigonometric S-matrices, especially in 
their RSOS form relevant in our case. 

We recall that our starting S-matrix has the form 

s ab (e) = x ab (e)s^(e)®s i ™\e). 

The total crossing-unitarizing prefactor Y ab is the product of CDD factor X ab and corresponding 
crossing-unitarizing prefactors er®, a^ for each S®(9) and S^ 1 . As was shown above, the 
correct TBA derivation, confirming the proper central charge in all the considered cases, was 
based on the assumption for the crossing-unitaring prefactor form of the full S-matrix (|17|) with 
the ip having the universal form for all the algebras (|42[). First, we start with a proof of this 
statement for A n and C n cases, when the full S-matrix sjh (6) is known for any pair a, b, 
since being free of multiplicities, the spectral decomposition, and crossing-unitarizing including 
prefactors ah , can be calculated explicitly. 

5.1 Vector- vector S-matrices and their fusion 

Recall the general RSOS structure of the S-matrix S^ b (6) for any algebra G n . In general one gets 
Sq)(9) as a result of fusion procedure starting from the fundamental S-matrix. By fundamental 
one means the S-matrix for fundamental representations, from which all other representations will 
appear after decomposition of tensor product of reps into irreducible ones. This representations 
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are usually called the defining representations. The defining reps for A n and C n algebras are 
their vector representations, while for B n and D n - their spinor ones. In the last case all other 
fundamental representations can be obtained as a tensor product of spinor ones. But, obtaining 
the vector representation in the tensor product of spinors, one can get from the vector represen- 
tation all the others by the same fusion procedure, except for the spinor ones. We don't know the 
explicit form of the spinor-spinor trigonometric RSOS S-matrices for B n and D n cases (one can 
guess they look quite complicated), although they are, as we said, the fundamental S-matrices 
in these cases. Instead, we will describe vector-vector S-matrices for all the cases. The struc- 
ture of RSOS S-matrix of the vector-vector representation is well known and described in the 
literature (see, for example, [24], [6]). We cite it here for completeness and reader convenience, 
following [25]. It is defined as scattering process of two kinks 

+ K cd (9 2 ) - K ab {6 2 ) + KuiOi) 

connecting different vacua a, b,c,d € A* of the theory from the weight lattice A* of the algebra 
G n . Weights of its vector representation are 



S = {ei -e ,...,e n+ i - e }, e = 

S = {0,±ei,...,±e n }, B n , 
£ = {±ei,...,±e n }, C n ,D n , 



S& 1 g 

n + l 



A. 



where ej is some orthonormal basis. Up to a prefactor scalar function Y this kink-kink S-matrix is 
proportional to the Boltsmann weights (BW) W of statistical lattice models with corresponding 
symmetry, constructed as solutions of Yang-Baxter interaction round the face equations [24] : 



0{) = Y{u)W 



d 



/(«; 



GbG c 



f(u)/2 



(59) 



where u = {62 — 6\)/m, f(u) = cu for some constant c, c — a,d — c,b — a,d — b G S, and G a will 
be defined below. The set of non zero BW for vector- vector representation looks as follows [24]: 



W 



W 



w 



a a + [i 
a + n a + 2/j, 

a a + \x 
a + fi a + /x + v 

a a + v 

a + (i a + fi + v 



sin(a; — Am) 



sin(aj) 



(60) 



sin a 



[IV 



+ An) 



sin(a^) 
sin(Au) / sm(a 



sin(a;) 



+ lo) sin(a 



[IV 



to 



sin 2 (a^) 



1/2 



{fi 7^ u) for A n case, and 
W 



w 
w 



a a + n 

a + fi a + 2(i 

a a + n 

a + /j, a + (i+ v 

a a + v 

a + fj, a + (i + v 



u 



u = 



u = 



sin(A — Ait) sin(a; — Ait) 

sin(w) sin(A) 
sin(A — An) sin(a M1 y + Ait) 

sin(A) sin(a^) 
sin(A — An) sin(Au) 



sin(a;) sin(A) 



sin(a^ + uo) sm(a^ u - lo) 
sin 2 (a^) 



1/2 



H ^ ±v 



(61) 
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w 



+s. 



fJ,U 



a a + v 
a + u a 

[A - An] [a^ u + u + An] 

[A] [Ctu-v + U] 



u 



sin(An) sin(a A1 _ i/ + uj — X + An) 



W 



a a 
a a 



sin(A) sin(a M _^ + a>) 
sin(A + An) sin(2A — An) sin(A — An) sin(An) 



sin(A) sin(2A) 



sin(A) sin(2A) 



Ja 



for B n , C n and D n , where u, v G X, a M = u(a + p) • a, (p is the sum of the fundamental weights 
of the algebra), ao = — w/2, = — a^, a^-j, = a M + a^, and the constants A are defined as 

Parameter i = 1 in A n ,B n , D n cases and t = 2 in C n case. We also used 



Ja 



a- 



_s{a ll + w) 



so,, 



n 



sin(a 



, v Q, G ao = i, 



/wk J 



Esin(a K +a;/2 - 2A) 
zr~r_ ; — — t*a K , 



sin(a K + uj/2) 



where a = —1 in the C n case and a = 1 in the B n ,D n cases. The function s(x) = sin(ix) in 
the B n ,C n cases, and s(x) = 1 for D n . The quantities G a used in are related to G afl by 
G afl = G a+fl /G a and can be written as 

n(+l) 

G a = e(a) ]^[ s(ai) sin(a; - a,-) sin(aj + a,-), 

i=l l<i<j<n(+l) 

where a = Yl'i^i^ a i e i defines a^, and e(a) is a sign factor chosen so that e(a + p)/e(a) = a. 
The models are called restricted since only the dominant weights 



a-6 <l 



(62) 



are allowed, where is the highest root of the algebra, and I is the level. 

The BW listed above satisfy a set of conditions important for the S-matrix construction: 



• unitarity 



where in A n case 



J2 W 



a e 
c d 



u W 



p(n) 



sin(a; — An) sin(a; + An) 
sin 2 (tj) 



(63) 



and in C n , D n cases 



p(n) 



sin(w — An) sin(u> + An) sin(A — An) sin(A + An) 
sin 2 (w) sin 2 (A) 



(64) 
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crossing symmetry (B n ,C n , D n cases): 



W 



a b 
c d 



W 



c a 
d b 



G b G c \ 1/2 
G a GdJ 



(65) 



crossing-unitarity relation (A n case) 



5> 



a e 
c d 



l + u)W 



a b 
e d 



1 - u 



. sin(A - An) sin(A + Ait) 

Obc ^Tr\ • ( 66 ) 

sin (w) 



The case A n is different, since vector representation is not conjugate to itself. 

Requirements of unitarity and crossing for the S-matrix (|59|). using (|63j) and Q66JI in the A n 
case, lead to the following functional constraint on the function Y(u), to which we will add two 
indices Y n i(u), explicitly emphasizing its dependence on rank n and level I: 



Yn,l(u)Y n7 i(-u) 
Y nj i{l + u)Y njl (l-u) 



sin 2 (w) 



s'm(u> + Ait) sm(uj — Xu) ' 

sm 2 (u;) 
sin(A + Ait) sin(A — An) 



The minimal solution of this system of functional relations (up to so called CDD umbiguities), 
which has no poles on the physical strip < u < 1, was found in [26]: 



Y n ,i{u) = exp<2 




in+l 



U 



x [I + n + 1] [n + 1] 



(67) 



, , n + 1 \ , / u(n + 1, 

(0 ( — u J - (n - 1 + 1) I — - - n - 1 



We recall that we use the short notations [z] = sinh(xz), (z) = cosh(xz). 

In the same way the functional relations on Y in the B n ,C n ,D n cases, using (|63 |l .(|65 |l . look 

like 



Y(u) 
Y(u)Y(-u) 



Y(l-u) 

p-\u), 



where p is defined in l|64j). The minimal solution of this system can be written in terms of Y n 



Y(u)=Y tg:tl (u)Y tgM (l-u) 



sin A 



sin a; 



(68) 



The S-matrices described here can be written in the spectral decomposed form with projectors 
onto irreducible representations appearing in the tensor product of two vector representations. 
(These projectors should be written in the IRF form). As was pointed out, the constructed 
vector-vector S-matrix is pole free and has no bound states. The bound states are produced by 
insertion of CDD factors, which have poles at the rapidity value corresponding to the desired 
projector in the spectral decomposition. It is well known that the closed and self consistent 
bootstrap procedure requires the pole in Su S-matrix in the channel corresponding to the second 
fundamental weight. Continuation of the bootstrap reproduces massive bound state multiplets 
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corresponding to all the fundamental representations of the algebra. This scenario was proved 
to be correct in the cases where the spectral decomposition of S a b is known for any a, b, and was 
conjectured to be correct in all the cases. One of the most effective methods of the construction of 
S-matrix spectral decomposition is the tensor product graph method (TPG). But unfortunately 
it works only when the spectral decomposition is multiplicity free. This is the situation in the 
A n ,C n cases where the tensor product of two representations with fundamental highest weights 
are multiplicity free. 

We briefly recall here the concept of TPG. First of all, there is requirement for the represen- 
tations, for which we build TPG, to be affinizable (for details see [27], [28]). All the fundamental 
representations of A n and C n are affinizable. Unfortunately, it is not the case for B n and Dri 
cases. TPG is a graph which is constructed by letting the irreducible components A and a of 
tensor product A a <g> A;, be the nodes of the graph joined by a link if A and a have of opposite 
parity and a C adj ® A. The parity of A is defined to be ±1 according to whether A appears 
symmetrically or anti-symmetrically in the tensor product in the limit q — > — 1. The concept of 
TPG works well when TPG is a tree and there are no multiplicities in the tensor product of two 
representations. In the Tables 1,2 one can find the TPG for any pair of fundamental weights of 
algebras A n and C n . 



Aq + Afe — » A a+ i + Afc_i ... — > A a+min ( TO+1 _ afe ) + Ab_ min ( n+1 _ ai b) 
Table 1. Tensor product graph for two fundamental representations a and b of A n (a>b). 



A a + At 


— ► A a+ i + A{,_i — > . 


• A a+ 5_x + Ai — > A a+ t 


I 


1 


1 


Aq-l + Ai_l 


— ► A a + \b-2 — ► • 


• A a+ 5_ 2 


I 


1 




Aq-6+1 + Al 






I 






A a -6 







Table 2. Tensor product graph for two fundamental representations a and b of C n (a > b). 
For a + b > n the graph truncates at the (n — a + l)th column. 

Given a TPG, the spectral decomposition of the R-matrix has the form 

Rab(x) = 

ft 

where the sum is over the irreps appearing in the tensor product decomposition, i.e. over the 
nodes of TPG, x is multiplicative spectral parameter, and the main rule dictated by TPG says: 
if there is an arrow from \i to v on the TPG then the coefficients and p v {x) satisfy 

p v {x) ~ a- V (/i)/2 - ' 

where 
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and the second Casimir /(//) = (fi + 2p) ■ fi. 

Using these TPG method one can construct the minimal S-matrix for scattering of any two 
fundamental multiplets of A n and C n [25]. It has the form: 

where 

^w=^wrirW»+ 2j " +M :;~ 6 ~ 2 ) (69) 

i=i*=i v n ) 

and 1" is defined by (|H7j) in A n case, and by (|H5|) - in C n case. Z ab coming from the R-matrix, is 
defined below. The R-matrices have the following form (6 > a) 

• A n case 

min(n+l— b,a) 

Rab(n) = Yl {-l) k+l p k ab (u)P Xb+k+K _ k 

k=0 

(by definition A n+ i = Ao = 0) with 



and Z a b in this case 



a 6—1 a 

^ab(«) = nn{2i + 2A;-a-6}[]{-2p-6 + a}. (70) 
j=i fc=i p=i 

Here and below we use a new notation 

sin(u;x/2 + Xu) 



{*} 



sin(oj) 



C„ case 



min(n— b,a) a—j 

jk _ {2p + b-a} -A- {2(n+l)+2q-b-a} 

PabW 11 {_ 2 p_6 + a } 11 {_ 2 (n + l) -2g + 6 + a}' 

P=l (7=1 



and Z„/, in this case 



(in A \ al ° a 
— ) TT + 2k ~ a ~ b}{a + b - 2(n + 1) - 2j - 2k} 

sin w j 
' j=i fe=i 



x [J{-2p - 6 + a}{-2(n + 1) - 2p + a + 6}. (71) 

p=l 
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• B„ and D„ cases 



In these cases affinizable among fundamental weights representations are only vector and 
spinor ones, and their tensor products are multiplicity free. We will consider here only S-matrix 
for two vector representations. It has the following spectral decomposition [23] [25]: 



S$(u) = {-2}{-g}Y gjl (u)Y g>l (l 



P2X 1 



{2} 
{-2} 



Px 2 + 



{2}{g} 
{-2}{-5} 



5.2 CDD factors 

The S-matrices presented above are crossing symmetric, unitary and minimal in the sense that 
they have no poles in the physical strip. As we said, in order to make bootstrap working, one 
should multiply these S-matrices by a set of properly chosen CDD factors. They appear as a 
result of fusion of the CDD factors for vector-vector S-matrices. It is useful to introduce the 
following notation for universal description of CDD factors 



X(a) = 




(We recall that t = 1 in A, B, D cases and t = 2 in C case.) 

• A n case 

Vector-vector S-matrix CDD factor is just A(2) and its fusion leads to the following CDD 
factor for a, b scattering 

a+b-1 

x ab ( U ) = n xu + i)x(j - 1). 

j=\a—b\+l, step 2 

One can see that there is a pole in X^u) at u = if a + b < n + 1, or at u = 2 — if 
a + b > n + 1. They correspond to particles a + b or a + b — n — 1 respectively in the direct channel. 
There is also a pole at u = corresponding to the particle \a — b\ in the cross channel. We 
will not discuss here the double poles, this discussion one can find in [25]. 

• C n case 

Here the vector- vector CDD factor has the form Xu = X(2)X(2g — 2) and its fusion gives 
a+b-l 

X ab (u)= H X(j + l)X(j-l)X(2(n + l)-j + l)X(2(n + l)-j-l). (72) 

j=|a— 6|+1, step 2 

As one can see from (|72jl. the pole structure here is more complicated - it exhibits poles up to 
the 4-th order. Discussion of the S-matrix analytic structure in this case one can find in [25]. 

• B n and D n cases. 

CDD factors for the vector- vector S-matrices in these cases has the same form An = X(2)X(2g— 

2). 
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5.3 Prefactor exponentialization 

Now we have to transform both X a t, and a^) to the exponential form in order to calculate 
the quantity Y a b = ^ (Xab&^ a abj use d in the TBA calculations. One can use for that the 

,2 



following identity 

sin(7ra) = exp 

valid for < a < 1, and also 



o 



dx2smh 2 (x (a - 1/2)) 
x sinh x 



sinh(A(g + ^a)) _ f S mh(^-2a)) \ 
sinh(A(0-i7ra)) " eXP \ x smh ^ x > sinh (§f ) j 

valid for < a < k. Straightforward but long calculations, using (jf)9^ . ((771|) . (|71j) . (|o"7|) . l|68j) lead to 
the following exponential representations valid in both A n and C n cases: 

X ab (u) = exp \2 J™ ^-e'™ (6 ab -2cothxA%(xj) S j , (73) 

*2(«) = exp ( 2 i~T e ~ ux1 ^ (74) 

where the kernels A^L are defined in Appendix and already appeared in the TBA calculations. 
It turns out that in B n and D n cases the vector-vector S-matrices fit the same general formulas 
(|73 |) .(|74 |) . One can see after some simple algebra, that these expressions give a universal answer 
for the quantity valid in all the cases described above 

where * means convolution. This coincides with the assumption made about the form of Y a b in 
the TBA calculations of the previous section. 



5.4 Special cases of S-matrices 

As we mentioned in the introduction, the form of S-matrix Q includes in it as subclasses S- 
matrices for other important two dimensional integrable models with Yangian symmetries. One 
of them are PCM (I, m — > oo). They are well studied, well defined for any G = A,B,C,D and 
their S-matrices are self consistent from the bootstrap point of view (see [23]). In particular, 
there are no double poles unexplainable by Coleman-Thun mechanism. The situation is more 
subtle with GN models (I — * oo,m = 1). Here the S-matrix conjecture 

S a b = X a t,Sa b , (75) 

wich was naively expected to be correct in all the cases G = A,B,D (see the footnote in the 
introduction about the C n case), was found to suffer from the "bootstrap violation" [23] in the B n 
case, while was shown to be correct in the A and D cases. Only recently the B n GN S-matrix was 
"corrected" [29] using elegant physical arguments about the symmetry of the Lagrangian. The 
presence of additional current (compared with D n case) makes the symmetry different, leading 
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to additional RSOS like vacuum degeneracy with additional kink structure. The answer for the 
S-matrix was shown to be 

Sab = X ab S ( a f ® S ab , (76) 

where S a b = 1 for any pair (a, b) except for (n, n), when S nn = Stci(^) - the RSOS S-matrix of 
the tricritical Ising model with rescaled rapidity. From the Lie algebraic point of view it is just 
Ai level 2 RSOS model described in the formulas (|60j). There are some identities and specific 
features of low level RSOS models. In particular, one can show that B n level 1 RSOS model is 
identical (up to a rescaling of the spectral parameter u) to A\ level 2 RSOS model. It is interesting 
to note that there are the same identities on the level of affine Lie algebras themselves, which 
were pointed out in the context of generalized parafermions in [31]. In this sense the S-matrix 
constructed in [29] naturally fits the general form of our S-matrix (fTj). The identity between ([75 j) 
and (UJ in A n ,D n cases follows from the fact that RSOS level 1 A n and D n S-matrices are equal 
to 1 and hence may be ignored in the tensor product of l|T|). 

The same situation one has with another class of integrable models - low level WZW models 
perturbed by current-current perturbation. In the most of cases they are equivalent to the GN 
models - just by fermionization of WZW models (see, e.g. [30]). They are well defined objects, 
while the GN model are not always well defined: for example, as we said C n GN model cannot 
be defined on the usual Majorana fermions, but level one C n WZW model with current-current 
perturbation is well defined and described by the S-matrix (^Q) with I — ► oo, m = 1. Here another 
interesting identity is valid. Analyzing Boltsmann weights together with the restriction condition 
(|62|). one can see the identity (up to a rescaling of the spectral parameter u) between C n level 
1 and j4 n _i level 2 RSOS models. This identity, which is absolutely analogous to the relation 
(B n )i~ (A\)2, was also pointed out in [31] for affine Lie algebras. 

Both of the identities express themselves also in the form of the TBA diagrams. Looking at 
the fig l.b. and at the fig 2. a, in the case m = 1, one can see that there remains one "spurious" 
node under the line of massive nodes on the fig l.b, and a chain of n — 1 "spurious" nodes under 
the massive line on the fig 2. a. They are nothing but the magnonic degrees of freedom, describing 
{A\)2 (TCI) model in the first case, and (A n _i)2 - in the second. This remark demonstrates how 
a correct TBA diagram could help to guess the correct S-matrix. 

In [25] it was shown, that S-matrix of the form l|75|) in the C n case sufferes from "bootstrap 
violation": there are double poles unexplainable by the Coleman-Thun mechanism. As we ex- 
plained, this confusion was related with a naive assumption about the S-matrix form (|75|). which 
was expected to be correct for current-current perturbation of the level 1 C n WZW model. The 
correct form of the S-matrix for this integrable model is S ab = X ab S^ ® and contains non 
trivial RSOS tensor factor. As we said, due to the isomorphism (C n )i~(^4 n _i)2 [31], this factor 
is physically natural. 

6 Discussion 

We derived TBA equations from the S-matrices (^Q) for all the infinite series if Lie algebras 
G = A n , B n ,C n , D n . We have shown that with assumption (|T7|) they give the Y-systems with 
the proper high temperature behavior reproducing the correct central charge of the cosets. The 
assumption Q17|) was shown to be correct in A n ,C n cases for any two fundamental multiplets, 
and also in B n ,D n cases for vector- vector multiplets scattering. 

It would be interesting to obtain crossing-unitarising prefactor for other B n , D n trigonometric 
RSOS S-matrices with available spectral decomposition - spinor-spinor, and vector-spinor ones, 
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in order to be sure in correctness of the assumption (|17[) also in these cases. 

Derivation of the TBA equations presented in this paper is a technical and quite straight- 
forward procedure. But we think it is not just necessary for completeness of the CFT - TBA 
relation picture. As we saw, one gets a feedback from this derivation procedure important for 
the form of the S-matrix itself. An attempt to reduce the derived TBA system to a known and 
studied Y-system, may give a hint for the correct form of S-matrix. For example, if one would 
start from the ([75 )) S-matrix for the B n invariant GN model, he will get the Y-system with one 
missing magnonic node on the TBA diagram, compared to the correct one (like Fig. lb). One 
immediately realizes what should be added in order to get the correct Y-system - tensoring of 
the S-matrix with RSOS S-matrix of TCI model will give a desired missing magnonic node. This 
simple logic may be useful in consideration of new, less studied integrable models. 
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8 Appendix 

8.1 Kernels and mass spectrum 



Here we list the matrices A G (inverse matrices for M a b/ [t } 
A n 

~A n [min(a, b)] [n + 1 — max(a, b)] 
[n + 1] 



~ An _ [imiiyu,, uji { ,v -r -l — UM^, U)\ 



~ R [min(a, b)] (n — ^ — max(a, b)) , v 

A% = [ - 1 h \ 2 ^Z, a ,6< n _i, (78) 

\B n <\B n _ [ Q ] n < „ _ 1 

2 (n - i) 

4(n-i)(l/2)- 



~ Cn _ [| min(a, b)} (| (n + 1 - max(a, 6))) 
(J (n + 1)) 



XD n _ [min(a, b)] (n - 1 - max(a, b)) 
(n-1) 



~D 1111111 L*. 17 / \ It J. 111CX-A-1 Uj. U ) . . 

Kb = ;J \„ r\ - J - LL , a,b<n-2, (80) 
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We also list here the kernels K Gn (their non zero elements), which are inverse for A Gn , and 
differ by a scalar factor from M^j [t~ b ] : 



K. 



ab 



Tab 



21 



a+l,b 



+ $a-l, 



K 



K b n , a,b<n-l, 
r Bn (1/2) 



(81) 
(82) 



nn—l 



K, 



nn 



A 
f 



n— In /j\ 

2(l/2) 2 

— (TT~ ' 



K 



ab 



f K ± (I) - «.!><n-L 



(83) 



K 



Cn 

nn—l 



A 



Cn 

n—ln 



(1/2) 
(1) ' 



K 



1. 



D r 



A 



A 



n—ln— 2 



K 



n— In— 1 



A^'\ a,fc<n-2, 
K 



(84) 



n—ln 
"■nn 



A nn-1 — - fv n-2n-l 



2(1)' 



1. 



As was mentioned above, mass spectrum vectors m a , with components corresponding to 
masses of different fundamental representation multiplets of algebras G n , form eigenvectors of 
matrices K Gn listed above with zero eigenvalues. Here we recall the mass spectra for different 
algebras 



A, 



m a = M sin ^ 



m, 



2M sin • 



■ a — a ^ 2n _ 1 

(a < n - 1) 
m n = M 



Msin 



2(n+l) 



D r , 



m n 



2Msin 

(a < n - 2) 

=m n = 



2n+2 

= M 



8.2 Solution of rf29j) . rt30 |) . 

First, we express p\~ x using the eq. (|2Tj) with j = 1: 

n-1 / 2Z-1 

Pi" 1 = E^n-l6^/2) 27iya 6 - £ A lfc (aV2)? 



.2(1) 
27 



~6 
A' 



fc=l 



Hence, using the identity 

a A Cn , A Cn 

At~x>m = f^ = 2(1/2)-^, 

we can write the last term in (|3U|) as 
Introducing notations 

na (V 2 ) / (0 a ~a\ ^ 1 / (0 n ~n\ n 

B = -^y (^i/ 2 CT ~ Px) ■> C = -^-y <7 -Pi)-D 
the system (|2T?|) . (|3U)) may be written as 



n-l 



6=1 1 ^ 

(V 2 ) + x n = a 



(1) 

Using the last equation, one can express x n and substitute it into the first one getting 

VlV- s 8 il/2) ^T b -B« + 6 (1/2) r 

The matrix in the l.h.s. has as an inverse the restriction of A F to the values a, b < n — 1, and 
we have the solution: 

^ > 6=1 

Using the obtained x n_1 in the second equation, we get 

n-l 

(1) ^ 

v ; 6=1 

where we used the identity 

4 c„ ( 1 / 2 ) 2 , i _ 4 c„ 

■"■n-ln-l Q\2 ' ^nn- 

Explicit use of the definitions of B a , C, D gives the expressions we used in the main text. 

8.3 Solution of (J3HD - flHED- 
We solve (f35|) with respect to x n 



T n _ A C„ r , (V 2 ) \ " A C n R 6 



2(1/2) W ^; / 



1/2 <J "l i f-^n^-l- 1 ) 
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and substitute it into We get an equation, which together with (J33J), can be written as 

E (*2* " <W V-iy (<n-i) J^ = fl- + ^-i^yC (85) 

where 

B " 2(1) I 1 
C = a^-p^-pl 

The inverse matrix for the one in the parenthesis in the l.h.s. of (|85|) coincides with A^ n , and we 
get 

n-1 

Z a = E^ b + ^fn-lC 
6=1 

One can now substitute x n ~~ l in the last form into x n and get finally the expressions (|39 |) (|4U |) . 
after collection of similar terms, using the kernel identity A®™_ x = 2(1/2)^4^, valid for a < n — 1. 
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